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ABSTRACT 

A new type of engine seal is being developed to meet the needs of advanced hypersonic 
engines. A seal braided of emerging high temperature ceramic fibers comprised of a sheath-core 
construction has been selected for study based on its low leakage rates. Flexible, low-leakage, 
high temperature seals are required to seal the movable engine panels of advanced ramjet- 
scramjet engines either preventing potentially dangerous leakage into backside engine cavities or 
limiting the purge coolant flow rates through the seals. To predict the leakage through these 
flexible, porous seal structures new analytical flow models are required. Two such models based 
on Kozeny-Carman equations are developed herein and are compared to experimental leakage 
measurements for simulated pressure and seal gap conditions. The models developed allow 
prediction of the gas leakage rate as a function of fiber diameter, fiber packing density, gas 
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properties, and pressure drop across the seal. The first model treats the seal as a homogeneous 
fiber bed. The second model divides the seal into two homogeneous fiber beds identified as the 
core and the sheath of the seal. Flow resistances of each of the main seal elements are combined 
to determine the total flow resistance. Comparisons between measured leakage rates and model 
predictions for seal structures covering a wide range of braid architectures show good agreement. 
Within the experimental range, the second model provides a prediction within 6 to 13 percent of 

the flow for many of the cases examined. Areas where future model refinements are required are 
identified. 
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NOMENCLATURE 

cross-area of seal 
yarn cross-sectional area 
fiber diameter 
gravitational constant 
mass flow rate of gas 
molecular weight of gas 
number of core yarns 
number of sheath yarns 
pressure downstream of seal 
pressure upstream of seal 
flow resistance of the seal 
Reynolds number 
universal gas constant 
absolute temperature 
seal dimensions (see Fig. 4) 
superficial gas velocity 
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y 0 half the clearance between the seal and its housing 

e porosity (see Eqs. (27) to (29)) 

<j> shape factor, defined in Eq. (4) 

6 braid angle (see Fig. 1) 

pf mass density of fiber 

p gas viscosity 

p gas density 

Subscripts: 
c core 

e edge 

si seal 

s sheath 

1,2. ..7 flow paths (see Fig. 4) 

INTRODUCTION 

Ramjet- scramjet engines require sliding panel seals to prevent combustion gases from 
leaking past the articulating engine panels, similar to articulating panel seals of turbojet two- 
dimensional converging-diverging nozzles. 1 However, new seals are required for advanced 
hypersonic engines because of higher thermal loads and the need to seal larger engine sidewall 
distortions. As a point of comparison, turbojet nozzle seals developed under the Augmented 
Deflector Exhaust Nozzle Program 2 used superalloy seals that sealed pressure differentials up to 
30 psi, sealed sidewall distortions up to 0.030 in., and were cooled to 1200 °F. Hypersonic 
engine seals, however, are required to operate at higher temperatures (1800 to 2000 °F), seal 
higher pressure differentials (up to 100 psi) and seal larger sidewall distortions (up to 0.150 in.), 
as described in Ref. 3. 
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A seal concept that shows promise of meeting these challenging demands is the braided 
ceramic rope seal being developed at NASA Lewis Research Center. The braided ceramic rope 
seal structure consists of a high density uniaxial core structure overbraided with an outer sheath 
for structural integrity, as shown in Fig. 1. Braided of emerging high temperature ceramic fibers 
this seal shows promise of operating hot and remaining flexible at temperatures up to 2000 °F. 
Early design studies identified important seal design parameters including: fiber diameter; yarn 
bundle size, fiber packing density, and percent core structure for low leakage. 

Accompanying the development of these engine seals, NASA is also developing engine seal 
flow models to predict the seal leakage through these porous seal structures. These seal flow 
models can be used during the design process in one of two ways: to predict performance losses 
associated with parasitic leakage through the seals; and to predict purge coolant flow rates 
through these seals where ambient engine flow temperatures exceed the seal’s operating 
temperature limit. The purpose of this paper is to provide an analytic means of predicting the 
gas flow through these braided structures and to determine quantitatively the relationship 
between gas leakage rate and the pore structure of the seal. 

THEORETICAL 
Definition of Flow Path 

As shown in Fig. 1, the flow path across a seal system can be divided into two categories: 

(l) flow through the seal and (2) flow around seal. The flow in the first path is related to the 
packing architecture of the seal itself and the flow via the second path is dependent on the 
surface properties of the seal and housing* 

E-foy through seal - Based on data from a large amount of experimental data obtained 
using a variety of packing materials, both spherical and granular in shape, Ergun 5 derived the 
following equation: 

~( P ° - P iK fcD) « 3 = 150(1 - .) j ?5 (1) 

p\l 2 t 1 - £ ( <j>D)pu/ft 


4 



where P denotes pressure, subscripts i the inlet and o the outlet. In the above equation, if 
the Reynolds number Re = (<£D) pu/fi (1 - e), is small (less than 10), then the constant 1.75 on 
the right side of the equation can be ignored. In the case of an engine seal, the flow is expected 
to be laminar as the gas leakage rate and the fiber diameter are small. This implies that the 
viscous term dominates in the above equation and the inertial term is negligible. Under such a 
condition, the pressure drop is proportional to flow velocity, u. 

In earlier studies the tortuous pore structure of the bed was modeled as a solid bed 
consisting of an assembly of capillaries with circular cross section. 6 The capillary model focused 
on the spaces or the pores in the porous solid. The best known equation proposed based on this 
approach is the Kozeny-Carman equation, which includes permeability coefficient as a function 
of porosity. One form of this equation is given below as: 

"( p 0 - Pi)*, 

u = ' ' 

150 

W e 3 

The shape factor, </> is defined as: 

^ __ area of sphere equivalent to particle volume 
actual surface area of particle 

The shape factor, <f> is unity for a sphere and 0.87 for a cylinder with its diameter equal to its 
length. The equivalent diameter of a particle is defined as the diameter of a sphere having the 
same volume as the particle. For a fiber with a diameter Dj and length L, the equivalent 
diameter is: 



D = 


1.5D? L 



( 3 ) 
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and the shape factor is: 


$ = (1.5) 2 / 3 


(L/D f ) 2/3 
L/D f + 0.5 


(4) 


Hence, (^D) can be expressed as: 


(^d) = i.5 n f L ^ Df 

L/D f -t- 0.5 


(5) 


If the ratio L/D f is very large the term (^D) will approach the value of 1.5D f . If direction of 
flow is across the axis of fiber, a situation that occurs in a seal containing substantial amount of 
longitudinal fibers, the length scale, L, in the above equation should be expected to be of the 
same order of magnitude as the diameter of the fiber. The parameter (<£D) can be thought of as 
characteristic dimension intrinsic to flow through the fibrous seal. 

The Kozeny-Carman equation predicts successfully the pressure drop in packed beds with 
porosity ranging from 0.3 to 0.6. For porous media with higher porosity such as most fiber beds 
and textile fabrics, a number of authors (e.g., see Ref. 8) have shown that predicted pressure 
drop is much greater than measured values. In the current application, for determining leakage 
rate of a gas through a seal having low porosity, Kozeny-Carman equation is a good starting 
point. Taking cross sectional area for gas flow as A f and the seal length as L, Eq. (2) can be 
rearranged as 


M _ ~( p ! - p ?) 

L snn'-y tL (■ - «)* 

A c e 3 (^D ) 2 

where ideal behavior of gas is assumed, and the density of gas p is based on an average value 
evaluated at the two end-point pressures as: 

, _ K + P >K , 

2R,T 
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Flow Around Seal .- Edge flow can be treated as a flow between parallel nonporous surfaces 


separated by a small gap. Assuming that the gap between the surfaces can be considered 
constant and equal to 2y o , one can relate the pressure difference across the seal to the gas 
leakage velocity as: 9 


3/iut 


( 8 ) 


Rearranging the above in the form of Eq. (6) gives: 


M 

r 


P‘ 


3MR g T t 

y 3 
J O 


(9) 


Flow Resistance .- Examination of Eq. (6) for flow though the seal and Eq. (9) for flow 
around the seal suggests the definition of flow resistance R as: 


R = 



( 10 ) 


M/L 

The flow resistance R is a function of properties of both fluid and seal architecture and, for this 
analysis, is assumed to be independent of the pressure difference across it and any compressive 
pressure the seal may be subjected to. 


Flow Modeling 

In a previous paper, 4 critical design parameters such as fiber packing density and fiber 
bundle size were identified through a combination of theoretical and experimental studies. Since 
the seal consists of both braids in the sheath and longitudinal fibers in the core, one expects 
different porosity values to be applicable in the two regions. In this investigation, two models 
are proposed to quantitatively evaluate the flow resistances. 

Model I .- As shown in Fig. 2, the seal is assumed to be a homogeneous fiber bed having a 
uniform constant porosity regardless of the core and sheath structures. Thus, only one value of 
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porosity is used to calculate the flow resistance. The gas leakage rate can be expressed as the 
sum of the leakages through the seal and around the seal, and is given by: 


m = K + Hi = ( p i - P o, 

L L L R 


where subscripts, si and e refer to the seal and edge, respectively. The individual leakages are 
given by: 


M. P i - P o 


Hi K - p2 


The flow resistances encountered in the flow path through the seal is determined from Eq. (6) 
and is given as: 

„ /x R_TtL (l — e) 2 

R„1 = 300 ? (14) 

Mwg c A c e 3 (0D) 2 

The edge flow consists of two parallel paths as shown in Fig. 2 and the flow resistance of each of 
these two paths may be summed in parallel as: 


R e = 


R el R e2 
R el + R e2 


R el = 9 - R,T 1 and R, 2 = 3^ill 
M.8. y 3 H.5, J 

where R e j is three times because of the longer path length, see Fig. 2. Since the edge 
flow and flow through the seal occur in parallel, the overall flow resistance of the seal system 
therefore given by: 


MR/r 


Re^sl 

R e + R 8 
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Model II - The second model, illustrated in Fig. 3, deals with a composite seal in which the 


sheath and core are allowed to have independent porosity values. The seal has a sheath and a 
core with porosities e s and s c , respectively. Flow resistances along the various flow paths 
illustrated in Fig. 4, are given as: 


R, = 9K__ 

3 


(17) 


R 2 = Rg = 300 K_ 


R 3 = R 5 = 300 


R 4 = 300 K 


(> - 4 

h 'K+dy 

k *2 t 1 ~ 

‘ •>»» 

(' ~ 


R 345 — R 3 + R 4 + Rq 


(18) 

(19) 

( 20 ) 
( 21 ) 


R 7 = 3K_ 


/iR.T 

where K = The flow resistance of the seal can be determined by summing the flow 

Mwgc 

resistances in parallel, given as: 


1 + 1 


R q Ro r 


+ JL 

Ri; 


(23) 


345 sx 6 


Flow resistance of edge flow is: 


R. 


RjRy 
Rj + R7 


(24) 


The total flow resistance of the seal for substitution into Eq. (11) is then given by: 

R e R s 

R = e 8 

R e + R» 


(25) 
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Calculation Basis 


Two important parameters in determining flow resistance through the seal are porosity and 
the characteristic dimension, D. The actual porosity of the seal in the application condition is 
expected to be lower than in the initial installed position. The porosity (method of 
determination discussed below) was assumed to be a constant for any given seal. With regard to 
D, since the bulk of the seal is made up of longitudinal fibers and that the number of fiber-fiber 
interfaces is significantly larger than the number of yarn-yarn interfaces, the characteristic 
dimension was taken as the fiber diameter, D^. The characteristic dimension (<^D), takes on 
values between 1.5D f and 0.75D f when the ratio L/D f is taken to be large or 0.5, 
respectively. The latter value gives better fit with experimental data under a variety of 
situations investigated. Therefore, all calculations presented in this paper are based on 
(^D) = 0.75D f . 

Another characteristic dimension is the distance, y Q in calculating seal edge leakage. In the 
present calculation, the clearance was assumed to be proportional to fiber diameter. Specifically 
y G is assumed to be 0.1 (^D). 

The cross section area of a yarn Ay (in. 2 ) can be determined from its denier and fiber 
density p { (g/cm 3 ) as: 


\ = 


yarn denier 
5.8xl0% f 


(26) 


where denier is a yarn density term used in the textile industry and is the yarn mass in grams 
per 9000 m of length. 

The following are the parameters used for calculation of various results reported: 


Pf = 2.54 g/cm 3 
R = 1.545 xlO 3 lb-ft/°R 
T = 528 °R 

g £ = 32.1 lb m ft/lb f sec 2 


D f = 10 pm (air) = 29 lb m /lb mole 

denier = 812 g per 9000 m (He) = 4 lb m /lb mole 

(0D) = 0.75D f p of air = 0.0175 cP 

y 0 = 0.1 (^D) p of helium = 0.019 cP 


10 


EXPERIMENTAL 


Braided Seal Specimens 

Eight seal specimens were made using 812 denier E-glass fibers (Owens Corning Glass, 
Granville, Ohio). The specimens were labeled Al through Hi and their architectural 
parameters, braiding angle, number of longitudinal yarns and number of braiding yarns, are 
summarized in Table 1. Specimens Gl and Hi have the highest number of longitudinal yarns 
while Al, Bl, and Cl have the lowest number of longitudinal yarns. 

Flow Measurement 

The experimental details of the flow measurement were described in an earlier paper by the 
authors. Seal specimens 1 ft in length were mounted in a specially developed test fixture and 
were leak tested under room temperature at various inlet pressure conditions in the range of 5 to 
80 psig. The pressure upstream of the seal was varied and the resulting leakage of gas (either 
air or helium) was measured. Lateral preloads were applied uniformly to the back of the seal 
with an inflatable rubber diaphragm at either 80 or 130 psig. The flow resistance of the seal was 
computed from the ratio of the difference of squares of absolute pressures over the mass leakage 
rate. 

Porosity 

An ultra-low viscosity embedding media (purchased from Polysciences Inc., Warrington, 

PA) was used as a rigidizing medium to infiltrate the specimen. Polymerization was 
accomplished at 70 C for 12 hr. The specimen then was cut and polished. Scanning electron 
micrographs were taken to determine the dimension of the seal cross section and packing 
geometry of fibers. The thickness of the braid sheath, t 2 , was measured using the electron 
photomicrographs. 
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Porosity of the seal for calculating flow resistances described earlier was obtained from the 
geometry of fiber layout and is given by: 

AJN C + N /cos (?) 

e = 1 - -J1 c 1 1 (2' 

t 2 

where N c and N s are the number of core and sheath yams and t 2 is the cross sectional area 
of the installed seal. Note that the seal is treated as a homogeneous fiber bed having a single 
average porosity value for Model I. 

The porosity of sheath and core sections of the specimens for Model II were determined 
from the following two equations: 


ff c = 1 


AyN c 


(28) 


6 * = 1 


AyN s /cos 0 
t 2 - 


(29) 


where t and t x are the overall width of the installed seal, and the width of the core region, 
respectively (see Fig. 4). 


RESULTS AND DISCUSSION 


Pressure Drop Correlations 

In Fig. 5, typical measured air leakage rates are plotted as a function of the difference of 
the squares of the pressure across the seal for specimens Al and Gl at preload pressures of 
80 and 130 psig. The linear relationship between the two variables is indicative of the validity 
of the pressure dependency presented earlier in Eq. (10). Although only two sample results are 
shown in Fig. 5, all eight specimens examined in this investigation showed excellent correlation 
with correlation coefficient lying in the range of 0.990 to 0.999. The slope of the line in Fig. 5 is 
equal to the inverse of flow resistance, 1 /R. 
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Flow Resistances for Different Gases 


Because of the many environments the seals are expected to operate in, it is important to 
be able to predict the flow resistance to various potential coolants or leakage gases. Shown in 
Fig. 6 is the measured resistance of helium plotted against the resistance of air for a wide range 
of seal architectures (specimens Al, Bl, Cl, Dl, Gl, and HI), pressure drop conditions (between 
5 to 80 psig), preload conditions (80 and 130 psig) investigated. If the seal’s pore structure is 
constant, flow resistance is directly proportional to viscosity and inversely proportional to 
density of the flowing gas (e.g., Eq. (6)). Hence, when we compare the flow resistance of helium 
to that of air in Fig. 6, we expect the slope of the straight line to be: 


SLOPE 


w 


Helium 




(30) 


The straight line indicated in Fig. 6 is the theoretical line with a slope of 7.9 obtained using 
Eq. (30). Experimental data given in Fig. 6 when fitted with a straight line yielded a slope of 
6.2, which compares favorably with the theoretical value of 7.9. 

Comparison of Measured and Predicted Leakage Rates 
The measured and predicted leakage rates for Al and Gl seals with widely different seal 
architectures are shown in Figs. 7 and 8 for applied pressure differentials up to 80 psi for both 
air and helium test gases. Also shown in the figure are the effects of lateral preload on seal 
leakage. Lateral preloads of 80 and 130 psig were applied to the back of the seal with a 
diaphragm compressing the seal against the adjacent sidewall. Key braiding and geometry 
parameters of these two seal structures denoted Al and Gl are listed in Table 1, 

Comparing the overall leakage rates between specimen Al and Gl one finds that the 
leakage rates for Gl are considerably less than Al. Specimen Gl meets the tentative leakage 
limit of 0.004 Ib/sec/ft (e.g., Ref. 3) for air pressure differentials up to 40 psi with a preload of 
80 psig. Specimen Al meets the leakage limit for pressure differentials only up to 30 psi. 
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In general, both models predict the leakage rates reasonably well with the measured data 
over the full pressure range for both air and helium test gases. Model II gives values closer to 
those measured at a preload of 80 psig, and Model I predicts values closer to those measured at 
130 psig. 

In examining Figs. 7(a) and 8(a), the versatility of Model II in predicting actual leakage 
rates is demonstrated. At a pressure differential of 40 psi and preload of 80 psig, the 
discrepancy between the measured and predicted leakage rates were only between 6 to 13 percent 
even though the overall leakage rates differed by a factor of 1.7. 

In all the cases, the fiber diameter was used as the basis of calculation. However, if yarn 
diameter is used as the equivalent diameter, the prediction of the gas leakage rate is very poor as 
it differs from experimental observation by more than four orders of magnitudes. 

Potential Sources of Modeling Discrepancy 

A potential source of the discrepancy between the measured and predicted leakage rates is 
the porosity dependence on preload. As the preload pressure is increased the fibers are urged 
closer to one another and closer to the adjacent surface making it more difficult for the air to 
flow around the fibers, thus increasing flow resistance. Neither of the models considered in this 
paper account for this porosity-load dependence and is presently under development. 

The choice of shape factor <j> has a considerable effect on the predicted leakage rates. For 
example, since in Eq. (6) the term containing ($D) is squared, increasing it from 0.75D f to 
0.80D f increases the predicted leakage by 14 percent. Selection of the shape factor <j> is based 
on the quality of the model fit to experimental observations and the seal porosity value used. 

As improved measured values of porosity become available, it is expected that the shape factor 
could change and approach the theoretical upper limit of 1.5. 
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SUMMARY AND CONCLUSIONS 


Two analytical models have been developed for predicting leakage rates of braided rope 
seals being developed for panels of advanced hypersonic engines. Both models are based on the 
Kozeny-Carman relations for flow through porous media, where the characteristic size dimension 
is a scaled fiber diameter (e.g., 0.75Df) based on experimental observations. 

The first model treats the seal as a homogeneous fiber bed having a single average value for 
its porosity. The second model treats the two-dimensional braided seal structures as a system of 
flow resistances analogous to a series of resistors in an electrical network. 

Based on the findings from the comparison between measured and predicted leakage rates, 
the following results were obtained: 

1. Leakage rates predicted using Model II agree favorably to the measured leakage rates for 
modest preloads (80 psig) for a wide range of braided seal architectures. Agreement within 6 to 
13 percent was observed at a pressure differential of 40 psi for seal specimens A1 and G1 whose 
overall leakage rates differed by a factor of 1.7. 

2. Theoretical predictions compared to experimental observations for air and helium 
indicate that relative resistance to leakage flow depends on the ratio of the quotients of each 
gas’s viscosity and molecular weight. 
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TABLE 1.— SEAL CONSTRUCTION DETAILS AND POROSITY DATA 


Sample 

number 

Braiding 

angle, 

or 

Core fiber, 
percent 

Average 

porosity a 

Thickness 
t d 

Core 

porosity* 5 

Sheath 

porosity c 

Al 

45 

39.6 

0.48 

0.11 

0.34±0.04 

0.54±0.02 

B1 

30 

39.4 

0.48 

0.12 

0.24±0.06 

0.56±0.01 

Cl 

10 

41.1 

0.50 

0.09 

0.50±0.03 

0.50±0.02 

Dl 

45 

57.2 

0.42 

0.09 

0.19±0.05 

0.58±0.02 

El e 

30 

60.7 

0.46 

n/a 

n/a 

n/a 

Fl® 

10 

55.0 

0.40 

n/a 

n/a 

n/a 

G1 

45 

82.2 

0.45 

0.02 

0.47±0.02 

0.37±0.16 

HI 

30 

81.3 

0.45 

0.02 

0.47±0.02 

0.32±0.18 


Calculated from Eq. (27). 

^Calculated from Eq. (28). 

Calculated from Eq. (29). 

d The accuracy of t 2 measurement was in the range of ±0.005. 
e El and FI were damaged during use. 
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Figure 1 .—Definitions of seal leakage flow paths. 
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Figure 2— Schematic diagram for 
flow Model I: uniform seal 
porosity. 


Figure 3. — Schematic diagram 
for flow Model II: independent 
sheath and core porosities. 
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Figure 4. — Flow paths used in calculating total flow resistance for 
Model li. 


Figure 5. — Measured seal leakage data validating seal leakage 
pressure dependence. 
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Figure 6. — Measured seal flow resistances of air and helium com* 
pared to theoretical expectations. 
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Figure 7. — Seal A1 ; measured leakage rates versus pressure drop 
(symbols) compared to predictions (lines). 



(a) Air. 



(b) Helium. 


Figure 8. — Seal G1 : measured leakage rates versus pressure drop 
(symbols) compared to predictions (lines). 
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